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Abstract 

We present the complete next-to-leading order (NLO) QCD corrections to tZ associated pro- 
duction induced by the model-independent tqg and tqZ flavor-changing neutral-current couplings 
Q_i. at hadron colliders, respectively. Our results show that, for the tuZ coupling the NLO QCD cor- 
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rections can enhance the total cross sections by about 60% and 42%, and for the tcZ coupling by 
about 51% and 43% at the Tevatron and LHC, respectively. The NLO corrections, for the tug cou- 



plings, can enhance the total cross sections by about 27%, and by about 42% for the teg coupling 



at the LHC. We also consider the mixing effects between the tqg and tqZ couplings for this process, 
which can either be large or small depending on the values of the anomalous couplings. Besides, 



the NLO corrections reduce the dependence of the total cross sections on the renormalization or 

O ■ 

factorization scale significantly, which lead to increased confidence on the theoretical predictions. 
And we also evaluate the NLO corrections to several important kinematic distributions. 
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I. INTRODUCTION 



The mass of the top quark is close to the electroweak(EW) symmetry breaking scale, and 
thus its decay and production at colliders are very important for the probe of the EW break- 
ing mechanism and new physics beyond the standard model (SM). Direct evidence of new 
physics at TeV scale may be not easy to find, while indirect evidence, such as modification 
of SM predictions originated from new physics interaction, is important as well. A good 
consideration is to investigate the single top quark production process via the anomalous 
flavor- changing neutral-current (FCNC) couplings. Within the SM, the FCNC couplings can 
only occur at the loop level, which are further suppressed by the GIM mechanism [lj. On 
the other hand, some new physics models (2,3] such as the two Higgs doublet models 4|, su- 
persymmetric models j^J, extra dimensions models js], may enhance these FCNC couplings 
to observable level. As the coming Large Hadron Collider (LHC) will produce abundant top 
quark events (about 10 8 per year), even in the initial low luminosity run (~ 10 fb"Vyear) 
8 x 10 6 top quark pairs and 3 x 10 6 single top quarks will be produced yearly, one may 
anticipate the discovery of the first hint of new physics by observing the FCNC couplings 
in the top quark sector. 

In general, any new physics at a high energy scale A can be described by an effective 
Lagrangian containing higher dimensional SM gauge invariant operators For the new 
physics induced top-quark FCNC couplings related to the gluon and Z boson, respectively, 
they can be incorporated into the dimension five effective operators as listed below {sj], 

9 z 

- 9s E -f^TVt 9 , + ih 9 tq ls)tGl v - — ^- J2 -f^ift* + <ls)tZ^ + H.c, (1) 

q=u,c ^ q=u,c 

where A is the new physics scale, T a are the Gell-Mann matrices, and Z^ v are the field 
strength tensors of the gluon and Z boson, respectively, kL (q = u, c) and rcf 9 (q = u, c) 
are real coefficients that define the strength of the couplings. And 9w is the weak-mixing 
angle, while / t q and h^ q are complex numbers satisfying |/ t q| 2 + |/^| 2 = 1 with V = Z, g 
and q — u,c. 

The CDF collaboration has set 95% confidence level (CL) limits on the branching ratio 
BR(t — > qZ) < 0.037 [9j, which corresponds to fc£/A < 0.908TeV _1 based on the theoretical 
predictions of t — > qZ at the NLO level in QCD 10h13|. The DO collaboration also provides 
a more stringent constrains, BR(t — > qZ) < 0.032 at a 95% confidence level pjj, which 
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corresponds to ftfj/A < 0.845TeV . Currently the most stringent experimental constraints 
for the tqg anomalous couplings are k? u /A < 0.013 TeV -1 and nf c /A < 0.057 TeV -1 , given 
by the DO collaboration [15JL and n 9 tu /A < 0.018 TeV" 1 and k 9 J A < 0.069 TeV" 1 given 
by the CDF collaboration [la], both based on the measurements of the FCNC single top 
production using the theoretical predictions, including the NLO QCD corrections 
and resummation effects 19], respectively. 

Since the observation of pp — > tZ is a clear signal of top flavor violation, and we do 
not know which type of new physics will be responsible for a future deviation from the SM 
predictions, it is necessary to study this process in a model-independent way. There are 
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already several literatures 20j, [2jJ discussing this process using the effective Lagrangian in 
Eq. (JTJ). However, they were either based on the LO calculations 20], or the NLO QCD 
effects are not completely calculated [2l|. So it is necessary to present a complete NLO QCD 
corrections to the above process, which is not only mandatory for matching the expected 
experimental accuracy at hadron colliders, but is also important for a consistent treatment 
of both the top quark production and decay via the FCNC couplings by experiments. In 
this paper, we present the complete NLO QCD corrections to tZ associated production via 
tqZ and tqg FCNC couplings with their mixing effects at hadron colliders. 

The arrangement of this paper is as follows. In Sec. [Ill we show the LO results for the 
process induced by tqZ FCNC couplings. In Sec. \\\\\ we present the details of the NLO 
calculations, including the virtual and real corrections. We discuss the process induced by 
tqg FCNC couplings and the mixing effects in Sec. HVl Sec. [V] contains the numerical results, 
and Section [VI] is a brief summary. 



II. LEADING ORDER RESULTS 



At hardron colliders, there is only one subprocess that contributes to the tZ associated 
production at the LO via the electroweak FCNC couplings, nf q : 

gq^tZ, (2) 

where q is either u quark or c quark. The corresponding Feynman diagrams are shown in 

Fig.m 

After sum over the spins and colors of the outgoing particles and average over the spins 
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FIG. 1: The LO Feynman diagrams for the single top quark production via the FCNC couplings 
without operator mixing. 

and colors of the incoming particles, the LO squared amplitude is 

32-K 2 aa.K 2 



\M B \ 2 (s,t) = -^^s^ (2m?-(3m^ + 4s + 2t)m t 6 + 

' ' 9q 3sin{2e w ) 2 K 2 s{t-m 2 f y 1 y z 1 1 

{2m 4 z - {2s + t)m 2 z + 2{s 2 + Us + t 2 ))m A t + {2m% - Atm\ 

+ {s 2 + 6ts + 5t 2 )m 2 z - 2t{3s 2 + 6ts + t 2 ))m 2 + t{-2m 6 z + 2 

{s + t)m z -{s + t) 2 m 2 z + Ast{s + t))) ) (3) 

where m t is the top quark mass, and m z is the Z boson mass, s, t, and u are the Mandelstam 
variables, which are defined as 

s = {pi+ p 2 ) 2 , t = (pi- p 3 ) 2 } u = {px- p A f. (4) 

After the phase space integration, the LO partonic cross sections are given by 

^ = U drWW ab - (5) 



2s . 

The LO total cross section at hadron colliders is obtained by convoluting the partonic cross 
section with the parton distribution functions (PDFs) Gi/p for the proton (antiproton) 

<y B = ^2 J dx ± dx z [ G a/p 1 {xi, fif)G b/ p 2 {x 2 , fi f )a^ b ] , (6) 

ab 

where ///is the factorization scale. 



III. NEXT-TO-LEADING ORDER QCD CORRECTIONS 

In this section, we present our calculations for the NLO QCD correcions to the tZ asso- 
ciated production via the electroweak FCNC couplings. The NLO corrections include both 
the virtual and the real corrections with the Feynman diagrams shown in Figs. [2J131 which are 



generated with FeynArts 



22j, and calculated with FormCalc 



23] . We use the dimensional 



regularization (DREG) scheme 24j with naive 75 prescription in n = 4 — 2e dimensions to 
regularize all the divergences. Moreover, for the real corrections, we use the two-cutoff phase 
space slicing method |25[ to separate the infrared(IR) divergences. 



A. Virtual corrections 
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FIG. 2: 1-loop Feynman diagrams for the single top quark production via the FCNC couplings 
without operator mixing. 



The virtual corrections contains both UV and IR divergences, with the UV divergences 
renormalized by introducing counterterms. For the external fields, we define all the renor- 
malization constants using the on-shell subtraction scheme 
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where C e = r(l + e)[(47r / u^)/m^] e and n/ = 5. For the renormalization of the strong coupling 
constant g S: and the FCNC couplings 5Z K a / A , we use the MS scheme [lo| : 

= ^r(l + e)(4vr)^ - J- + ^LC7 e — , 
gs 4vr 1 n ; V 3 2 7 e uv 12vr 6 Cl/v ' 

rv 1 

SZ Kfg/A = ^r(l + e)(4vrr — , (8) 
and the running of the FCNC couplings are given by [10] 

A A W^)/ 



(9) 



with /3o = 11 -2n//3. 

All the UV divergences cancel each other, leaving the remaining IR divergences and the 
finite terms. Because of the limited space, we do not shown the lengthy explicit expressions 
of the virtual corrections here. The IR divergence of the virtual corrections to the partonic 



total cross section can be factorized as 
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b7r ^ e IR \ m t m t m t * J e IR J 

where -D e = [(47r/i^)/s] e /r(l — e). In order to cancel these divergences, we need to extract 
the IR divergences in the real corrections, which will be shown in the following subsection. 



B. Real corrections 

The real corrections consist of the radiations of an additional gluon g q — > t Z g, or 
massless quark(anti) in the final states, g g — > t q Z, q q(q, q') — > t q(q, q') Z, q'q' — > tqZ 
as shown in Fig. |3j It should be noted that in our NLO calculations of the process induced 
by tqZ couplings, we did not include the contributions from the SM on-shell production of 
the top quark pair with subsequent rare decay of one top quark, pp(p) —>ti—>t + q + Z, and 
also the corresponding interference terms, following the diagram removal scheme proposed 
in reference This procedure does violate the gauge invariance because certain diagrams 
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We have also 



are removed, but the influence to the numerical results is small as show in 
crosschecked the numerical results by using another method preserving gauge invariance [29] 
where an invariant mass cut of the Z boson and light quark is adopted. And the results of 
the invariant mass cut method agree well with the ones of the diagram removal method if 
we require the invariant mass to be out of the range of ±17r t of the top quark mass, where 
T t is the width of top quark. 




FIG. 3: Feynman diagrams of the real corrections for the single top quark production via the 
FCNC couplings without operator mixing. 

1. Real gluon emission 

For real gluon emission, the phase space integration contains both soft and collinear 
singularities. We adopt the two-cutoff phase space slicing method to isolate all the IR 



singularities [25| , which introduces two small parameters 5 S and 5 C to divide the phase space 
into three parts. The soft cutoff 5 S separates the phase space into the soft region E 5 < 5 s y/s/2 
and the hard region, 

a R = a R + a s , (II) 
furthermore, the hard piece can be divided into two sub-regions by S c , 



a 



H_~HC + ~HC (12) 



The hard noncollinear part <r HC is finite and the phase space integration can be calculated 
numerically. For the soft region, in the limit that the energy of the emitted gluon becomes 
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small, i.e. E 5 < 5 s y/s/2, the amplitude squared ^2\M(qg) — » tZ + g)\ 2 can be factorized 
into the Born amplitude squared times an eikonal factor $ ci k 



J2\M(qg)^tZ + g)\ 2 ^ (47ra s ^)^|M s | 2 $ cik , 
where the eikonal factor is given by 

Ca s 1 m 2 — u 

2 (pi ■ Ps)(P2 ■ Ps) 2C A (Pi ■ P5)(Pi ■ Pn) 



(13) 



; ik 



Ca m 2 — t 



m; 



(14) 



2 (p 2 -Ps)(P4 'PS) (P4-P5) 2 ' 

where C4 = 3,CV = |. Moreover, the three-body phase space in the soft limit can also be 
factorized 



dT 3 (qg ^tZ + g)^% dY 2 (qg -4 tZjdS. 
Here dS is the integration over the phase space of the soft gluon which is given by 
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The parton level cross section in the soft region can be expressed as 
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After the integration over the soft gluon phase space[25j, Eq. (fl7|) becomes 
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where A and B±, are given in Appendix. 
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In the hard coUinear region, E$ > 5 s \fs/2 and — 5 c s < t i5 < 0, the emitted hard gluon is 



coUinear to one of the incoming partons. As a consequence of the factorization theorem 
3"i| the matrix element squared for qg — > tZ + g can be factorized into the product of the 
Born amplitude squared and the Altarelli-Parisi splitting function 

■2P qq {z,e 



J2\M(qg^tZ + g] 



1 2 collincar 



{A7ia s ^)J2\M B \ 2 



+ 



(20) 



where z denotes the fraction of the momentum of the incoming parton carried by q(g) with 
the emitted gluon taking a fraction (1 — z), and the unregulated Altarelli-Parisi splitting 
functions are written explicitly as [25 ] 

■1 + z 2 



P qq (z,e) = Cv(y 



2N(- 



— z 

Z 1 — 2 



+ z(l-z) 



(21) 



V 1 — z z 

Moreover, the three-body phase space can also be factorized in the coUinear limit, for 
example, in the limit — 5 c s < £15 < it has the following form (25] 

(4tt) 



dTs(qg —ttZ + g) co ii^|. ar dT 2 {qg — > tZ; s' = zs 



-dzdt 15 [-(l - z)t 15 ] 



(22) 



Thus, after convo 
region is given by 



167r 2 r(l - e) 

uting with the PDFs, the three-body cross section in the hard coUinear 



25] 
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+Pgg(z,e)G g / p (x 1 /z)G q / p (x 2 ) + (xi 4+ x 2 )] 

where G q ( g y p (x) is the bare PDF. 



-)* c ' [P q q(z,e)G q/p (x 1 /z)Gg /p (x 2 ) 
dz ( \ — z s 



dx\dx 2 . 



(23) 



2. Massless (anti)quark emission 

In addition to the real gluon emission, a second set of real emission corrections to the 
inclusive cross section for pp — > tZ at NLO involves the processes with an additional massless 
q(q) in the final state. Since the contributions from real massless q(q) emission contain initial 



state coUinear singularities we need to use the two cutoff phase space slicing method [25( to 



isolate these collinear divergences. The cross sections for the processes with an additional 
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massless q(q) in the final state can be expressed as 
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where 



^ \ a u {qa -)■ tZ + q{q))G q/p (x 1 )G a / p (x 2 ) + 

(a=q,q,q') 
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(24) 



P gg (z,e) = P ffff (z) = -[z 2 + (1 - z) 2 \ - z(l - z)e, 
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Pgq(z,e) - Pgq{z) - C F [- - _ 
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(25) 



The a c terms in Eq. (l2~4"j) represents the noncollinear cross sections for the q(q, q, q') and gg 
initiated processes which can be written in the form 



o C = y s { £ \M(q(q, q, q')) <™T tZ + (q, q, q')\ 2 + \M(gg) tZ + g) | 2 } ^,(26) 

(a=q,q,q') 

where dT^ is the three-body phase space in the noncollinear region. The other terms in Eq. 
(124")) are the collinear singular cross sections. 



3. Mass factorization 



After adding the renormalized virtual corrections and the real corrections, the parton level 
cross sections still contain collinear divergences which can be absorbed into a redefinition 



of the PDFs at the NLO, namely through mass factorization 



32 



331 ] . This procedure, in 
practice, means that first we convolute the partonic cross section with the bare PDF G a / P (x) 
and then use the renormalized PDF G a / P (x, /i/) to replace G a / p (x). In the MS convention 
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the scale-dependent PDF G a / P (x, jif) is given by [25 ] 



G a / P (x,n f ) = G a / P (x) + ^2 \~~ ) 

P ^ e ' 



a s T(l - e) 
27rrfl - 2e) 



x 



^4 



" x dz 

x / —P a p(z)Gp/ p (x/z). 



(27) 



Then the 0(a s ) expression for the remaining collinear contribution can be written in the 
following form: 
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{G q /p{Xl, flf)Gg/p[X2, ///) + G q / p (Xl, {J,f)Gg/p{X 2 , flf) 



+(xi -H- x 2 )}rfxi<ia;2, 



G q / p {x 1 ^ f )Gq/ p (x 2 ^f) 



(28) 



where 



with 



^ c (<? = C F (21n5 s + 3/2), 

Af(</ -> </</) = 2JVln < y, + (lliV-27i / )/6, 

.-1— 5 S <5 

G a/p (x,nf) = Yl 
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dy 

y 



G^/p{x/y,fi f )P a p{y), 



P a ,{y) = P a p{y) in(5 c — 4r) - P^(y), 



(29) 
(30) 
(31) 
(32) 



(33) 



y 

where N — 3,n/ = 5. 

Putting together all pieces of the real correction, we can see that the IR divergences from 
the real correction can be written as 



R . a s „ r 13 /„, . s , , .m? — w. , .m? — i, 43 \ 1 1 o . . 

* = + 1 4111 ^' + H ^r ] ~ 91n( ^r> + t ) ^K- (34) 

and all the IR divergences from the virtual corrections in Eq. (TlOT) are canceled exactly, as 
we expected. 
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IV. CONTRIBUTIONS FROM THE ELECTROWEAK AND STRONG FCNC 
COUPLINGS AND THE MIXING EFFECTS 



couplings, nf . However, for the tZ associated production process, there are additional 
contributions from the strong FCNC couplings, kL, and the mixing effects between these two 



In previous sections, we only consider the contributions from the electroweak FCNC 

z 

v;\rvi + In o L' f iv wi iv HI W I ^mmlmiTD L ■ 

tq ) 

couplings. Since the magnitudes of the coefficients k^XV = Z,g) depend on the underlying 
new physics, the mixing effects may be significant in certain model. The 0(a s ) corrections 
to the process pp — > tZ induced by tqg are similar to ones induced by tqZ, so we don't 
show its analytical results, and only present the mixing effects in this section. 

At the LO, the contributing Feynman diagrams are show in Fig. HJ and the squared 
amplitudes are present in the Appendix. The NLO corrections, which include the loop 
diagrams and the real emission diagrams, are shown in Fig. [12] and Fig. 






FIG. 4: The LO Feynman diagrams for the single top quark production via the FCNC couplings 
with operator mixing. 



The relevant renormalization constants are the same as ones in Eq. ([7]) and ([8]), except 



that we introduce additional renormalization constants. We adopt the definition in Ref. 
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£,.„■ + <>4,r =( -,>-)( 1 + 5Z -^ Z || | (30) 

<5Z Zg 1 + 6Z Z z 

where the operators {i = g,Z) are defined as O g = gsQ cr ' J ' u T a (f^ q + ihfqj^tG", Oz = 
sde^Q alXU (ft<i + ihtq^tZ^, and 5Z gg = 5Z K 3 q/A , 5Z ZZ = 8Z R zj K . At the 0{a s ) level, 
SZ K z / A is presented in Eq. ([8]), and other renormalization constants are given by: 
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«,z = ^r(i + £ )(4.).^, 

on e uv 
6Z Zg = 0, (36) 

SZgz is defined at the level of the cross section, and ci, C2 are defined as follows: 

s 3 — Qfsin 2 6w 



sinO 



w 



c i = g*zL9g L Qf sin °w - g Z R9 gR - 

n o * s 3 -Q f sin 2 e w 
c 2 = 9zR9 gR QfSin6 w - g ZL g gL — — , 

where Qf is the electric charge of the fermion, and S3 is its third component of the SU(2)l 
gauge group. 



The renormalization group running for k v are modified to 



-V 



3/3 , 



A A 



+ _*£^iZ( —sin 2 8 w -4) (77^0 -J7 3 ft>) , (37) 



A A ' A 

where 77 = a " ) . 

The IR divergence appearing in the virtual corrections has the same form as given in 
section IHIt except using the LO amplitude including the contributions from both the elec- 
troweak and the strong FCNC couplings instead of a B '. 

V. NUMERICAL RESULTS 

A. Process via tqZ FCNC couplings without operator mixing effects 

We first consider the tZ associated production via the tqZ FCNC couplings, including the 
NLO QCD effects on the total cross sections, the scale dependence, and several important 
distributions at both the Tevatron and LHC. All the SM input parameters are taken to 
be |34|: 

m t = 172.0GeV, a s (M z ) = 0.118, a = 1/128.921. (38) 
And we set the electroweak FCNC couplings, allowed by current experiment, as follows: 

K Z Jk = K Z Jk = 0.5TeV-\ (39) 
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The running QCD coupling constant is evaluated at the three-loop order [35J and the 



CTEQ6M PDF set [36( is used throughout the calculations of the NLO (LO) cross sec- 
tions. Both the renormalization and factorization scales are fixed to the sum of the top 
quark and the Z boson mass. 

In Table HJ we list some typical numerical results of the LO and NLO total cross sections 
for the tZ associated production via the electroweak FCNC couplings. It can be seen that, 
for the tuZ coupling the NLO corrections can enhance the total cross sections by about 
60% and 42%, and for the tcZ coupling by about 51% and 43% at the Tevatron and LHC, 
respectively. 



FCNC coupling tuZ (LO) tuZ (NLO) tcZ (LO) tcZ (NLO) 
LHC (pb) 15.9 22.5 1.29 1.85 

Tevatron (fb) 55.5 88.6 1.62 2.45 

TABLE I: The LO and NLO total cross sections for tZ associated production via the electroweak 
FCNC couplings at both the LHC and Tevatron. Here /i = rat + m 2 , K tq/-^ = 0.5TeV _1 . 

In Figj5] we show that it is reasonable to use the two cutoff phase space slicing method 
in our NLO QCD calculations; i.e., the dependence of the NLO QCD predictions on the 
arbitrary cutoffs 5 S and 5 C is indeed very weak. While the Born cross sections and the 
virtual corrections are cutoff independent, both the soft and collinear contributions and the 
noncollinear contributions depend strongly on the cutoffs. However, the cutoff dependence 
in the two contributions (a s + o coU and a HC + a c ) nearly cancel each other, and the final 
results for a NLO are almost independent of the cutoffs. We will take 5 S = 10 -4 in the 
numerical calculations below. Generally S c being 50 — 100 times smaller than 5 S is sufficient 



for accurate calculations to a few percent |25j, so we take 5 C = 5 s /50 in our calculations. 

In Figs. |6]we show the scale dependence of the LO and NLO total cross section for three 
cases: (1) the renormalization scale dependence \i r = /i, fif = m t + mz] (2) the factorization 
scale dependence /i r = m t + m Z) fif = fi; and (3) total scale dependence \i r = = //. It 
can be seen that the NLO corrections reduce the scale dependence for all three cases, which 
make the theoretical predictions more reliable. Fig. [7] give the pt distributions of the top 
quark and the Z boson, respectively, and Fig |8] shows the invariant mass distributions of 
the Z boson and the top quark. 
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FIG. 5: Inclusive total cross sections for pp — > tZ + X at the LHC as a function of 5 S in the phase 
space slicing treatment. The 5 C is chosen to be 5 C = 5 s /50. 




FIG. 6: Scale dependence of the total cross sections for the gu initial state subprocess at both 
LHC and Tevatron, the black lines represent the LO results, while the red ones represent the NLO 
results. 



15 




200 300 
PT.,„( GeV > 




200 300 

p TZ (GeV) 



FIG. 7: pt distributions of the top quark and the Z boson at LHC, the black and red lines represent 
the LO and the NLO results of the FCNC single top production, respectively. 




FIG. 8: Invariant mass distributions of the Z boson and the top quark, the black and red lines 
represent the LO and the NLO results of the FCNC single top production, respectively. 



B. Process include operator mixing effects 

In this subsection, we present the numerical results of the tZ associated production 
via the electroweak and strong FCNC couplings, including the NLO QCD effects and the 
mixing effects. We investigate the NLO QCD effects on the total cross sections and the scale 
dependence at the LHC. For the numerical calculations, we take the same SM parameters 
as above subsection. 
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In Table [TT1 we list some typical numerical results of the LO and NLO total cross sec- 
tions for the tZ associated production via the strong FCNC couplings, assuming «f q /A = 
O.OlTeV -1 , allowed by current experiment. It can be seen that, the NLO corrections can 
enhance the total cross sections by about 27% for the tug couplings, and by about 42% for 
the teg coupling at the LHC Nevertheless, the contributions to the total cross sections from 
the process induced by the tqZ FCNC couplings are still dominate. 

FCNC coupling tug (LO) tug (NLO) teg (LO) teg (NLO) 
LHC (fb) 141 180 7.6 10.8 

TABLE II: The LO and the NLO total cross sections for the tZ associated production via the strong 
FCNC couplings at the LHC. Here u = m t + m z , «£ q /A = O.OlTeV -1 , and \g g L? = \g g R? = L 

After considering the mixing effects, the total cross sections of the tZ associated produc- 
tion via FCNC couplings can be factorized as: 



- = ^ gL | 2 4 f + Mm\ 2 ^ff + A,Cf ) 2 + 

[AJt*(9*M - A 5 Re(g sK g* ZR )] (40) 
where A4 represent the contributions from different couplings and mixing effects. And, their 
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FCNC coupling tuV (LO) tuV (NLO) tcV (LO) tcV (NLO) 
LHC (fb) 147 188 8.1 11.5 

TABLE III: The LO and the NLO total cross sections for the tZ associated production via the 
FCNC couplings at the LHC. Here fi = m t + m z , k^/A = O.OlTeV -1 and \g g L? = \g g R\ 2 = 
Re(9gLg* ZL ) = Re(g gR g* ZR ) = 1 . 

numerical expressions at the LHC can be written as 

= {ri5|^ L | 2 (^) 2 + 699|^ R | 2 (^) 2 + 63.7(^) 2 

9 Z 

+ [7Q.0Re(g gL g ZL ) - 8A.lRe(g gR g ZR )]^^}(ph • TeV 2 ), 

(41) 

a**, = {925|, gL | 2 (^) 2 + 874|, gR | 2 (^) 2 + 90(^) 2 

+ [W7Re(g gL g ZL ) - lURe(g gK g* ZR )] ^}(pb ■ TeV 2 ), 

(42) 

Jll = {38.9| %L | 2 (^) 2 + 37.2|^ R | 2 (^) 2 + 5.15(^) 2 

+ [6.97Re(g gL g ZL ) - 7.47ife(^ R )] ^}(pb ■ TeV 2 ), 

(43) 

<lo = {56.7| %L | 2 (^) 2 + 51.5|^ R | 2 (^) 2 + 7.38(^) 2 

9 Z 

+ [9.13Re(g gL g ZL ) - 9.83i2e(^ R )] ^}(pb • TeV 2 ), 

(44) 

where gn,giR are chiral parameters: 

9n = /tq - ih\ v g iR = / t q + i/ijq, \g iL \ 2 + |^iii| 2 = 2. 

In Table [TTll we list some typical numerical results of the LO and NLO total cross sections 
by choosing a special set of parameters (for simplicity, we set fi^/A = ^ q /A = O.OlTeV -1 ) 
and fix n = m t + mz- For the g u — >■ t Z , the NLO corrections can enhance the total cross 
sections by about 28%, and for the g c — >■ t Z process, by about 42% at the LHC. 

To investigate the contributions from the operator mixing effects, we present the counter 
curves for the variables Kz/A, /%/A, and Re(g* L gzL), Rz{.9gr9zr) in Figs. [HI [TO] 
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FIG. 9: The contour curves of the total cross sections versus the parameters K g /A and k z /A for 
the process induced by tuV and tcV FCNC couplings. Here we set \g g L\ 2 = |%r| 2 = R e {9gL9*zh) = 
Re(g gR g* ZR ) = 1. 
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FIG. 10: The contour curves of the total cross sections versus the parameters Re{g g L x gzh) and 
R- e (ggR x 9zr) for the process induced by tuV and tcV FCNC couplings. Here we set \g g L,\ 2 = 
\g gR \ 2 = \gzh\ 2 = \gzn\ 2 = 1 and Kg /A = k z /A = O.OlTeV- 1 . 
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FIG. 11: Scale dependence of the total cross sections for the gu initial state subprocess at the 
LHC, the black lines represent the LO results, while the red ones represent the NLO results. 

In Figs. [TT] we show the scale dependence of the LO and NLO total cross section at LHC 
for three cases: (1) the renormalization scale dependence n r — /i, /i/ = m t + mz] (2) the 
factorization scale dependence fi r = m t + m z , /i/ = /i; and (3) total scale dependence \i r = 
fif = fi. It can be seen that the NLO corrections reduce the scale dependence significantly 
for all three cases, which make the theoretical predictions more reliable. 

VI. CONCLUSIONS 

We have calculated the NLO QCD corrections to the tZ associated production via the 
tqZ and tqg FCNC couplings at hadron colliders, respectively, and we also consider the 
mixing effects of these two couplings. Our results show that, for the tuZ coupling the 
NLO corrections can enhance the total cross sections by about 60% and 42%, and for the 
tcZ coupling by about 51% and 43% at the Tevatron and LHC, respectively. The NLO 
corrections can enhance the total cross sections by about 27% and 42% for the tug and 
the teg couplings, respectively, at the LHC. The mixing effects between the tqZ and tqg 
FCNC couplings for this process can be either large or small depending on the values of 
various anomalous couplings. If we set n g /A = k z /A = O.OlTeV -1 and \g g i\ 2 = \ggn\ 2 = 
Re(g g Lg* ZL ) = Re(g g Rg ZR ) = 1, the NLO corrections can enhance the total cross sections 
by about 28% for tuV couplings, and by 42% for tcV couplings at the LHC. Moreover, the 
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NLO corrections reduce the dependence of the total cross sections on the renormalization 
or factorization scale significantly. 
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APPENDIX 

In this appendix we give the LO results of the tZ associated production induced by the 
tqg and tqZ FCNC couplings, and the definition of the A and B±. 
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= w^mn-s - - <M m ^ {s - m ' )2(2m - - (3 

m 2 + 4s + 2t)m^ + (2m 4 - (2s + t)m 2 + 2(s 2 + Ats + t 2 ))m 4 + (2 

m^ - 4tm 4 + (s 2 + 6ts + 5t 2 )m 2 - 2t(3s 2 + Qts + t 2 ))m 2 t + i 

(-2m* + 2(s + t)m\ - (s + t) 2 m 2 + 4st(s + t)))K 2 z m 2 z + 9t 2 

<7zr(s - m t 2 ) 2 (2mf - (3m 2 + 4s + 2t)m 6 t + (2m 4 z - (2 

s + t)m 2 z + 2(s 2 + 4ts + t 2 ))m 4 + (2m* - 4tm 4 + (s 2 + 6ts + 5t 2 ) 

m 2 - 2t(3s 2 + 6ts + t 2 ))m 2 + t(-2m 6 z + 2(s + t)m 4 - (s + t) 2 m 2 + 4st 

(s + t)))K 2 z m 2 z + 6st# gL #zL(s - m 2 )(t - m 2 )(-2(4sw 2 - 3)(3t - m 2 z )m 6 t + 2(s 

(3 - 4sw 2 )m 2 + t(s(16sw 2 - 9) + 4(4sw 2 - 3)t))m 4 + t(-8s 2 sw 2 + 2(3 - 4sw 2 ) 

t 2 + (-8tsw 2 + s(16sw 2 - 9) + Qt)m 2 z + 12s(l - 2sw 2 )t)m 2 + st 2 (8ssw 2 + 

8tsw 2 + (3 — 8sw 2 )m 2 — Qt))K g K z m 2 + Qstg g ^gz^{s — m 2 )(t — m 2 ) 

(8sw 2 (m 2 — 3t)m t 6 + (2t(16tsw 2 + s(16sw 2 - 3)) - 8ssw 2 m 2 )m 4 - t 

((8sw 2 - 6)s 2 + 6(4sw 2 - l)ts + 8sw 2 t 2 + (8tsw 2 + s(3 - 16 

sw 2 ))m 2 )m 2 + st 2 (8tsw 2 + (3 - 8sw 2 )m 2 + s(8sw 2 - Q)))K g K z m 2 z - s(t - m 2 ) 2 

((2(3 - 4sw 2 ) 2 m 2 (m 2 - t)m\ + 2(3 - 4sw 2 ) 2 m 2 (t(s + 2t) - 2(s + t)m 2 z )m\ - 

2(3 - 4sw 2 ) 2 (s 2 + Ats + t 2 )m 2 z (t - m\)m\ + st(-4(3 - 4sw 2 ) 2 

(s + t)m 4 + (2s 2 (3 - 4sw 2 ) 2 + 6t 2 (3 - 4sw 2 ) 2 + 3s(64 

sw 4 - 64sw 2 + 15)t)m 2 - 9st 2 )m 2 + s 2 t 2 (2(32 

sw 4 - 24sw 2 + 9)m 4 - 2(32sw 4 - 24sw 2 + 9)(s + t)m 2 z + 9s 

t))9& + 3g R (32sw 4 m 2 (m 2 - t)m* + 32 

sw 4 m 2 (t(s + 2t) - 2(s + t)m 2 z )m\ - 32sw 4 (s 2 + 4ts + t 2 ) 

m2 z(t ~ m l) m t + st(-64sw 4 (s + t)m\ + (32s 2 sw 4 + 96t 2 

sw 4 + 3s(64sw 4 - 32sw 2 + 3)t)m 2 - 9st 2 )m 2 + s 2 t 2 

(2(32sw 4 - 24sw 2 + 9)m 4 - 2(32sw 4 - 24sw 2 + 9) 

(s + t)m 2 z + 9st)))K 2 g )}, (45) 
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81n -g±i 

A = y h , (46) 

P+ 0T / j (gij - + g) + m f 2 + - g - 2m \ ( ^Ij-J^T . 
L " 2Ll2( 2^T^ J+2Li2 _ -t 



In 2 ( '? ' ) ~ ^ [ ) , (47) 



2 (m t 2 - u) / 2 V - 1 



B = -2L12 f.-(^lK+^-lK-^ + » + ^ +2Li2 / + 
V 2{ml-s-u) ) 2 V 



where /3 = t /l - / 2 4sm 2 ^ n 2 - 
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FIG. 13: Feynman diagrams of the real corrections for the single top quark production via the 
FCNC couplings with operator mixing. 
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